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ABSTRACT 



IMPULSE RADIO ULTRA- WIDEBAND 
COMMUNICATION OVER FREE-SPACE OPTICAL 

LINKS 



A composite impulse radio ultra-wideband (IR-UWB) communication system is 
presented. The proposed system model aims to transmit UWB pulses over several kilo- 
meters through free-space optical (FSO) links and depending on the link design, the 
electrical estimates of the FSO system can be directly used or distributed to end-user 
through radio-frequency (RF) links over short ranges. However, inhomogeneities on 
the FSO transmission path cause random fluctuations in the received signal intensity 
and these effects induced by atmospheric turbulence closely effect the system perfor- 
mance. Several distinct probability distributions based on experimental measurements 
are used to characterize FSO channels and using these probabilistic models, detection 
error probability analysis of the proposed system for different link designs are carried 
out under weak, moderate and strong turbulence conditions. The results of the anal- 
ysis show that depending on the atmospheric conditions, system performance of the 
composite link can have high error floors due to the false estimates of FSO link. The 
system performance can be improved by employing error control coding techniques. 
One simple solution employing a convolutional encoder and Viterbi decoder pair is 
also analyzed in this thesis. Another important system parameter that is the average 
channel capacity of the FSO system is analyzed under weak and moderate turbulence 
conditions. Theoretical derivations that are verified via simulation results indicate a 
reliable high data rate communication system that is effective in long distances. 
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OZET 

OPTiK SERBEST-UZAY KANALLARINDA 
DARBE-RADYO ULTRA-GENi§BANT HABERLE§iMi 

Bu gah§mada karma darbe-radyo ultra-geni§bant haberlesme sistemi 6nerilmi§tir. 
Onerilen sistem modeli ultra-geni§bant sinyalleri optik serbest uzay (OSU) iizerinden 
uzun mesafelerde ta§imasi amaglanmi§tir ve sistem tasarimma bagh olarak elektriksel 
OSU kestirimleri direkt olarak kullamlabilir veya son kullamciya radyo frekansi (RF) 
iizerinden kisa mesafelerde dagitim yapabilir. Ancak, OSU gonderim kanallarmdaki 
tiirde§sizlikler alman sinyalin yogunlugunda dalgalanmalara neden olur ve atmosfer- 
deki tiirbiilans kaynakh bu etkiler sistem performansim yakmdan etkilemektedir. OSU 
kanallari tammlamak igin deneysel olgiimlerden gikarilan belirli istatistiksel dagihmlar 
mevcuttur ve bu modelleri kullanarak onerilen sistemin hata olasihk analizi zayif, orta 
ve sert tiirbiilans bolgelerinde incelenmi§tir. Yapilan teorik analiz sonuglarma gore 
atmosfer ko§ullarma bagh olarak sistem performansmda yanks OSU kestirimlerinin ge- 
tirdigi yiiksek hata tabanlari gozlemlenmi§tir. Sistem performansmm iyile§tirilmesi igin 
hata kontrol kodlamasi 6nerilmi§tir. Ornek olarak, OSU sisteminde gonderici tarafmda 
evrisjimli kodyahci ve ahci tarafmda Viterbi kodgoziicusii kullamldigi durum incelenmi§ 
ve bu durumun hata olasihk analizi yapilmi§tir. Onemli bir ba§ka sistem parametresi, 
ortalama kanal kapasitesi OSU kanallar igin zayif ve orta tiirbiilans ko§ullarmda in- 
celemi§tir. Simulasyon sonuglari ile dogrulanan teorik gikarimlar, onerilen sistemin 
uzun mesafelerde giivenilir bir §ekilde yiiksek veri hizlarma gikabilecegini gostermi§tir. 
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1. INTRODUCTION 



Ultra-wideband (UWB) systems characterize those that employ very narrow pulses 
in time usually on the order of nanoseconds to convey information. Consequently, these 
pulses occupy very large bandwidths in frequency domain. Due to their characteristics, 
UWB systems have unique attractive features introducing new advances in several wire- 
less communications fields such as radar, imaging, networking and positioning systems. 
These aforementioned fields benefit from several advantages offered by UWB systems 
such as enhanced ability to penetrate through obstacles such as concrete walls, ultra 
high precision ranging, possibility to transmit at very high data rates and support 
multi-user, relatively small size and less power consumption [TJ. Until 2001, the use 
of UWB systems were limited solely to military applications and with the regulations 
introduced by the Federal Communications Commission (FCC) in the United States, 
these new features of UWB systems have gained its attraction in several areas [2j. The 
FCC in the United States has allocated a huge unlicensed frequency spectrum ranging 
from 3.1 GHz to 10.6 GHz for UWB systems. These pulses need occupy at least 500 
MHz of bandwidth within the allowed spectrum or a fractional bandwidth of more than 
20 percent where the fractional bandwidth defines the ratio of bandwidth B over the 
center frequency fc- Since UWB systems do not employ any carrier to transmit the 
signals, the definition of center frequency needs to be clarified. The lower and higher 
frequencies of —10 dB emission points on the power spectrum of the transmitted signal 
are denoted by fi and /# and their arithmetic mean determines the center frequency 
of the signal, fc = (/l + /h)/2. Also, the bandwidth of the signal is defined as the 
difference of these two frequencies, B = f H — f L [T]. 

Although wireless propagation channels have been studied for more than 50 years 
and several mathematical models are introduced in literature to describe their propa- 
gation characteristics, the FCC in the United States have accepted the IEEE 802.15.3a 
channel model for the indoor and outdoor UWB channels with the contribution from 
Intel Corporation in 2002. This regulation for UWB systems are primarily developed 
for Wireless Personal Area Networks (WPAN) and to enable the coexistence of UWB 
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systems with the current technologies, namely, the Bluetooth and the IEEE 802.11 
wireless local are networks (WLAN), the FCC in the United States has defined the 
spectral mask that every UWB device needs to accompany which is presented in the 
next chapter. 

Although UWB systems bring their unique advantages, they are limited to short 
distances, usually up to 10 meters, due to severe attenuation of these narrow pulses. To 
overcome this short range limitation in radio frequency (RE) links, free-space optical 
(FSO) systems have emerged as a cost efficient solution. Compared to RF systems, 
FSO links offer much larger bandwidth and capacity, less power consumption, more 
robust against eavesdropping and better protection against interference [3]. Although 
FSO systems are mainly used in inter-satellite and deep space communications, these 
systems have attracted considerable attention to solve the last-mile solution in the past 
decade [3|. 

Contrary to other emerging alternatives to RF links, theoretical background of 
FSO systems are known and they are very similar to the fiber-optical communication 
systems. In these systems, information can be conveyed in the intensity, frequency, 
phase or polarization of the signal. Among these attributes, intensity is typically used 
to carry information due to its simplicity and this type of modulation employing varia- 
tions in intensity is called intensity modulation. Because of the complexities associated 
with phase or frequency of the received signal, these type of modulation are not pre- 
ferred in most FSO systems [5]. 

Along with the power regulations concerning eye safety, FSO system performance 
are also limited by the atmospherical effects. The inhomogeneities in the atmosphere 
cause fluctuations in the amplitude and phase of the received signal over long dis- 
tances and these effects severely degrade the performance of the system. Besides the 
inhomogeneities in the atmosphere, aerosol scatterers such as rain, snow and fog and 
building-sway as a result of wind loads, thermal expansion and weak earthquakes are 
the other factors degrading the system performance EJ- For instance, under severe 
fog conditions, it has been reported that the signals attenuate on the order of hundreds 
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of decibels per kilometer [7]. 

1.1. Research Overview and Contributions 

This work addresses the short range coverage problem of UWB systems that op- 
erate in RF links. To overcome this bottleneck, an alternative technique is proposed 
and UWB signals are transmitted through FSO links over several kilometers in the last 
mile. In FSO systems, the impairments caused by the atmospheric turbulence closely 
effect the system performance on the transmission path and therefore, there has been 
extensive studies and models proposed to describe these conditions and these can be 
found in [S]- |T3] - Recent studies using these models investigate the performance of 
FSO links under atmospheric turbulence effects such as those in [U [5j El [HI [15] . In 
their pioneering work, Zhu and Kahn have applied maximum-likelihood sequence de- 
tection (MLSD) technique to mitigate the effects of atmospheric turbulence at the cost 
of increased complexity They have also included maximum-likelihood detection 
under the case when spatial diversity is available. In [6], the authors have analyzed 
the performance of FSO system for both correlated and independent channels under 
weak turbulence conditions. The performance analysis of a FSO system that use pulse 
position modulation scheme under weak turbulence conditions is presented in [H] and 
in their system model, the authors employ avalanche photodiode to receive the signals. 
In his study, Kiasaleh carries out the same analysis for the same type of photodetector 
but this time, the analysis includes the effects of moderate turbulence conditions [I]. 
Another recently published work in [15], the authors investigate the system perfor- 
mance under a more accurate statistical distribution to model the moderate weather 
conditions. Besides error analysis, the ergodic channel capacity of FSO links in the 
weak turbulence regime is derived in [16] and for the moderate regime, its closed-form 
expression is presented in [IT] . Also, Hranilovic et. al derive outage capacity analysis 
of an FSO system considering the effects of pointing errors in [IS] and Uysal et. al 
apply a relay system to their FSO system model in [19j . 

However, this work focuses solely on conveying UWB signals through FSO links 
over long distances and to this end, several system models are proposed. In the first 
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system model, the generated optical UWB signals can be transmitted over FSO links 
over long distances and at the receiver the detected electrical outputs can be used 
right away or distributed to the end-user within the building using ethernet or fiber 
cables at very high-speeds. It should be emphasized again that FSO systems can 
achieve very high data speeds that are compatible with fiber optical and ethernet 
cables can support. As an extension of the first case, the second scenario considers 
the case that the UWB signals are again conveyed over FSO links in the last-mile 
but this time the electrical outputs are distributed to the end-users by the RF link 
in the last-step by employing a simple UWB transmitter and receiver pair. Third 
scenario uses error-control techniques in the FSO link to reduce the possible errors 
induced by the turbulence effects. Although under good weather conditions the system 
performance can achieve superior results, as the link visibility decreases or weather 
conditions gets worse, error-control coding is a necessity to endure the availability 
of the link and improve the quality of service of the system. For all three scenarios, 
corresponding system models are presented and their analytical results for the detection 
error probability under different turbulence conditions are analyzed. These theoretical 
results are verified via simulations that confirm the presented results. To sum up, the 
biggest contribution of this thesis is the complete error analysis of the aforementioned 
system designs that enable long distance transmission of UWB signals through FSO 
links. This contribution is valuable because UWB systems suffer from their short 
range limitations and the proposed models in this thesis would help to overcome this 
bottleneck. 



1.2. Organization of the Thesis 



Chapter [2] provides brief information about impulse-radio ultra-wideband sys- 
tems. It presents several factors that are considered to model channel distributions in 
RF links and introduces the spectral mask proposed by several regulatory institutions 
in Europe and United States that enables the coexistence of these devices with the 
other equipments in the same frequency spectrum. In its first section, several sources 
of fading are discussed in detail and particular examples are given from daily life. In the 
next section, the accurate distribution to model the RF-UWB links are discussed. This 
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accepted model is mainly based on Saleh-Valenzuela model with some modifications. 
Parameters definitions and suggested values of these parameters are also presented 
towards the end of the chapter. Chapter [3] focuses on wireless optical systems. It dis- 
cusses several commonly used modulation techniques in these systems and introduces 
some conditions on the reception of the transmitted signal along with some important 
parameters that determine the link quality. Its first section explains the characteristics 
of FSO channel in detail and the last section of the chapter introduces the statistical 
models for different turbulence conditions. Chapter @] discusses the proposed system 
model in this thesis. Along with the devices used in the system model, their mathemat- 
ical interpretations are also given. Detection error probability analysis of the several 
possible scenarios carried out. For instance, closed form expression are derived for the 
scenario where the FSO link is used to convey the information over several kilometers 
and the RF link to deliver the signal to the end-user in the last-step and simulations 
are made to verify theoretical results. This chapter is concluded with the channel ca- 
pacity analysis of the FSO links. Finally, the last chapter addresses several concluding 
remarks and suggestions for possible future work. 
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2. IMPULSE-RADIO ULTRA- WIDEBAND SYSTEMS 



In wired channels such as coaxial cables, wave propagation is predictable and usu- 
ally stationary whereas in radio channels, the behavior of wave propagation extremely 
random in nature [20] . There are several sources that effect electromagnetic wave prop- 
agation in radio frequency (RF) channels such as absorption, diffraction, reflection and 
scattering which all gives rise to multiple paths. On the path between transmitter 
and receiver, the objects are highly likely to be absorb, diffract, reflect or scatter the 
energy of the waves causing the signals to travel by various amplitudes, phases and 
paths before it reaches to the receiver. The propagation of these waves (echoes of the 
original transmitted signal) are called multipath propagation and the phenomenon of 
amplitude and phase fluctuations in the received signal due to multipath propagation 
is called multipath fading [21] . Since each wave travels different lengths of paths, pulses 
arrive at the receiver at different times with different amplitudes and phases. There- 
fore, to describe the behavior of the channel in different environments, several channel 
models are proposed. These proposed models are mostly based results measured from 
experimental data. 

Several multipath components arrive to the receiver and collectively determine 
the amplitude and phase of the received signal. The first statistical model which is 
still widely used is Rayleigh fading channel model and it applies to the conditions 
when the bandwidth of the signal is small and the delays of the individual multipath 
components do not interfere with each other in the received signal model. This signal 
model assumes that there are almost no significant scatterers nor dominant signal 
reflectors in the channel. The other widely accepted channel model is Ricean fading 
channel and it considers the case when the channel exhibits fixed scatters or signal 
reflectors in the channel [22]. 

Although these models are sufficient to characterize channels for narrow-band 
systems, more accurate description to capture the behavior of the multipath com- 
ponents is needed to characterize systems that occupy ultra-wide bandwidths. For 
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this purpose, several studies have been made. In 2002, the IEEE has established a 
standardization group, IEEE 802.15.3a to a develop standard for UWB Personal Area 
Networks (PANs). They were to select the multiple access schemes, modulations and 
a commonly agreed channel model along with the other necessity of the standard. To 
this end, IEEE 802.15.3a Task Group had been formed and established a standard 
channel model to be used for the evaluation of PAN physical layer proposals. The 
Federal Communications Commission (FCC) in the United States has also decided on 
the data rates for the standard model under the regulated power levels below —41.3 
dBm in Part 15 rules to coexist with the already available services such as Global 
Positioning System (GPS) and the IEEE 802.11 wireless local are networks (WLANs). 
IEEE 802.15.3a standard model aims to transmit at data rates up to 100 Mb/s at 10 
meters, 200 Mb/s at 4 meters and higher data rates at smaller distances |23j . 

The allocated spectrum to transmit these signals ranges from 3.1 GHz to 10.6 
GHz and this ultra-wide bandwidth allows for several bandwidth-demanding low power 
applications in wireless communications such as short-range high-speed Internet access, 
vehicular radar and sensors, asset and personnel tracking, imaging through steel and 
walls, surveillance and medical monitoring applications pp. It should be emphasized 
that UWB signals are on the order of nanoseconds scale so that their power spectral 
density occupies such bandwidths without the need for any carrier. Depending on 
the pulse width, these signals can be designed to occupy all the 7.5 GHz spectrum 
in a single-band or several multi-band UWB signals can be created such that each 
occupies at minimum bandwidth of 500 MHz. The former is called as impulse-radio 
ultra-wideband (IR-UWB) signalling and the latter is referred to as multi-band UWB 
signalling where each has its own consequences and this work focuses solely on IR-UWB 
signalling. 

On the other hand, UWB technology faces several challenges such as bit error 
rate, capacity and throughput limitations due to low-power regulations, network flexi- 
bility, immunity to synchronization and many more. Along with these challenges, the 
FCC has released a spectral mask that limits the equivalent isotropic radiated power 
(EIRP) spectrum density of the UWB pulses to fully coexist with the other aforemen- 
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tioned technologies and this is shown in Fig. 12.11 and Tables 12.11 and 12.21 along with 
the regulations introduced by European Conference of Postal and Telecommunications 
(CEPT) and International Telecommunications Union (ITU). Pulse shaping filters are 
commonly used to generate such pulses and several optimal pulse shapers satisfying 
the UWB spectral mask are presented in [24|. 



Table 2.1. CEPT regulations for indoor and outdoor UWB equipments 





/ < 3.1 (GHz) 


3.1< / < 10.6 (GHz) 


/ > 10.6 (GHz) 


Indoor (dBm) 


-51.3 + 87 log (//3.1) 


-41.3 


-51.3 + 87 log (10.6//) 


Outdoor (dBm) 


-61.3 + 871og(//3.1) 


-41.3 


-61.3 + 87 log (10.6//) 



Table 2.2. ITU regulations for indoor and outdoor UWB equipments 



Frequency (GHz) 


Indoor (dBm) 


Outdoor (dBm) 


0.96 - 1.61 


-75.3 


-75.3 


1.61 - 1.99 


-53.3 


-63.3 


1.99-3.1 


-51.3 


-61.3 


3.1 - 10.6 


-41.3 


-41.3 


Above 10.6 


-51.3 


-61.3 




5 10 15 5 10 15 

Frequency (Hz) x 10 s Frequency (Hz) x 10 s 

Figure 2.1. Several indoor and outdoor spectral masks in which UWB equipments 

need to satisfy. 
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In the following sections, several sources of multipath fading in radio frequency 
(RF) channels and IEEE 802.15.3a standard channel model is be discussed in detail. 

2.1. Sources of Multipath Fading 

2.1.1. Reflection and Propagation 




Figure 2.2. Reflection of a wave incident on a medium. 

This phenomena happens at an interface between two different (dislike) media. 
When waves pass from one medium to another, the wave front suddenly changes its 
direction from which it originated. This is called reflection. Reflections may be specular 
(mirror-like) or diffuse (i.e. retaining only the energy, not the image) depending on the 
nature of the interface. For instance, when waves pass from a dielectric medium to a 
conductor medium (dielectric-conductor), the phase of the reflected wave is retained 
whereas depending on the incident angle, it may or may not be retained when waves 
pass from a dielectric medium to a dielectric medium (dielectric-dielectric). 

2.1.2. Diffraction 

The other phenomena that propagating waves experience is diffraction. It is 
the spreading out of waves. All waves tend to spread out at the edges when they 
pass through a narrow gap (for distances that are comparable to the wavelength of 
the electromagnetic wave) or encounters an object. Instead of saying that the wave 
spreads out or bends round a corner, it is commonly referred as the wave is diffracted 
around the corner. 
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Intuitively, the longer the wavelength of a wave, the more it will diffract. Another 
important remark is that the diffraction loss increases with increasing frequency. A 
familiar example from daily life is the reflection of image on a compact disc (CD) 
or digital versatile disc (DVD). Here, the surface of a CD/DVD acts as a diffraction 
grating for the reflection of images incident on its surface. This is the main reason that 
when one looks at the back of a CD/DVD, the tracks of it act as a diffraction grating 
and produce iridescent reflections of the sunlight. Another example from daily life is 
the color of spider webs. 

Specular 




Figure 2.3. Diffraction of a wave from a medium with curvature ends. 
2.1.3. Scattering 

The obstacles or inhomogeneities in the medium force to deviate the propagation 
of light from its straight trajectory. When the light is scattered by one large scat- 
terer, then it is called single scattering and similarly, when there are several sources 
of scattering in the medium that interact with the light, then it is named as multiple 
scattering. Mathematically, single scattering is treated as random since the location of 
the single scatterer in the propagation path is not known whereas multiple scattering 
is considered to be deterministic since the random nature of several scatterers can be 
averaged out. An example for single scattering is when an electron is sent to an atomic 
nucleus. The exact position of the nucleus relative to the electron is not known and 
this bring the random nature of this interaction. To solve these type of problems, in 
most cases, probability distributions are used to describe single scattering. A light 
passing through a thick fog can be given as an example for multiple scattering. Several 
models are proposed in literature to model scattering coefficients depending on the en- 
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vironment [25] such as Rayleigh scattering where the inhomogeneities in the medium 
are small compared to the wavelength of the light, Mie scattering where both sizes 
are comparable to each other and non-selective scattering where the obstacle size is 
much larger than the wavelength of the light. Rayleigh scattering is the reason that 
the sky is blue in the day time and reddish in the sunset. Mie scattering occurs on the 
lower portions of the atmosphere where the clouds overcast the sky and the abundant 
particles are large. Non-selective scattering is the reason that the fog and clouds look 
white since each color of the light is scattered around in equal portions resulting in 
white color. 




■ 



Figure 2.4. Random scatterers in a medium effecting the propagation path of light. 
2.2. The IEEE 802.15.3a Standard Channel Model 

This model aims to capture the multipath characteristics of typical indoor en- 
vironments where IEEE 802.15.3a devices operate. In addition to the UWB channel 
measurements performed mainly in 2002, a number of measurement campaigns were 
carried out by the participants of the task group [26j. Throughout these campaigns, 
three channel models were considered in general: Tap-delay line Rayleigh fading model, 
the Saleh-Valenzuela (S-V) model and the A-K model and these are summarized in 
[23] and [27]. 

The experimental measurements indicate that the transmitted signals through 
multipath RF environments arrive to the receiver as clusters. The proposed UWB 
channel model is based on the cluster approach first proposed by Turin and others 
in 1972 [28] and further formalized by Saleh and Valenzuela (S-V model) which was 
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proposed in 1987 [29] only with one significant modification. It was experimentally 
observed that when lognormal distribution is used to model multipath gain magnitude 
instead of Rayleigh distribution, it would better fit to measurement data. 

The S-V model which assumes the multi-path components of the same pulse 
arrive as clusters, the time of arrival of these clusters are modelled as Poisson arrival 
process with rate A 

pC^-i) = Ae- A &- T *- l) (2.1) 

where Tj and are the arrival times of the j th and (J — l) th clusters, respectively. 
The arrival time of the first cluster is usually set to zero, T\ = 0. Similarly, within each 
cluster, consecutive multi-paths arrive with a Poisson arrival process with rate Ai 

p{r jk \r u -i) k ) = Aie- Al ^- T «-^ (2.2) 

where and 77j-i)fc are the arrival times of the j th and (j — l) th paths (rays) of the 
k th cluster, respectively. As again, the arrival times of the first rays within each cluster 
is set to zero, Tj\ = for j = 1,...,N. Having defined the arrival rate processes, 
the proposed IEEE 802.15.3a multipath model has the following discrete-time impulse 
response 

N K(j) 

h(t) = K ^2^2 a jkS(t-Tj-T jk ) (2.3) 

j=l k=l 

where k denotes the log-normal fading random variable, and K(j) denotes the 
number of observed clusters and multipaths within each cluster, respectively. Tj is the 
delay of the j th cluster and {oijk} denotes the channel coefficient of the j th cluster k th 
multipath. The channel coefficients can be further can be represent by 



Oijk — PjkPjk 



(2.4) 
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where pjk is a Bernoulli random variable taking equiprobable values of {±1} and it 
takes account for the signal inversion due to reflection. The other random variable, (3jk 
is the log-normal distributed channel coefficient that belongs to the j th cluster of k th 
path. This coefficient is expressed as 

P jk = 1 K W 20 (2.5) 

where Gaussian random variable with mean fijf. and variance cx 2 fc and 

further represented as 

Kjk = [ijk + £jk + Cjk (2.6) 

where £jk and Qk are two Gaussian random variables representing the fluctuations of 
the channel coefficient on each cluster and on each contribution, respectively. The 
variances of these two random variables are represented by <r| and cr 2 , respectively. 
The first term on the right hand side in (12.61) . fijk is closely related to the exponential 
power decay factors of amplitude of the clusters and multipath by 

E [\P jk \ 2 ] = E [|lO^+^+^)/ 20 | 2 ] = E [lAx,! 2 ] exp _ T ±\ (2 . 7 ) 

^ln( J E;[|/3 00 | 2 ])-(^ + ^)] lnl (af + a 2 ) 



10 

Hjk = — 



In 10 20 



(2A 



where _E[|/3oo| 2 ] is referred to as normalization factor of the total received power and 
sometimes denoted by Qq for simplicity. Often in communications, in order not to 
amplify or attenuate the transmitted signal, the channel coefficients are normalized to 
unity for each realization in the system model and it requires that 

N K(j) 

= 1 ( 2 - 9 ) 

j=l k=l 



Once again, let us remark that the arrival times of each cluster and each multipath 
within a cluster is modelled by two Poisson process which were previously defined in 
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(ED) and (EZD - 

The amplitude gain k is also assumed to be a log-normal r. v. with the relation 
IO 9 / 20 where g is a normal random variable with mean g and variance a 2 g . The 
mean g depends on the average total multi-path gain G and expressed as 

lOlnG _<&*W 
y0 In 10 20 V 1 

where ln(-) stands for the natural logarithm. G is dependent on average attenuation 
exponent ip by 

G = G /D^ (2.11) 

where Go is the reference power gain evaluated at D — 1. A = 10 \og w (ETx / E^x) 
is the path loss at a reference distance Dq = 1 in dB. In literature, this is sometimes 
denoted as PL and it is related to reference value for power gain G by 

G = 10- Ao/1 °. (2.12) 

In 2003, based on several experimental results, Ghassemzadeh and Tarokh suggested 
values for Aq = 47 dB and ip = 1.7 to be used in line-of-sight (LOS) environment 
and Aq = 51 dB and ip = 3.5 in NLOS environment [30]. Using these definitions, the 
parameters for IEEE 802.15.3a UWB channel model can be summarized in Table [2731 

Since different in-door environments exhibit different channel characteristics and 
the channel parameters need to be updated. These different characteristics arise from 
several facts. Depending on the distance between the transmitter and receiver and 
whether a LOS component exists or not several channel model scenarios are been 
proposed. The first scenario, channel model 1 (CM1) is valid for distances of — 4 
meters with LOS component. The second scenario, CM2 is valid for the same distances 
but it considers NLOS cases. For distances of 4 — 10 meters with NLOS cases, CM3 
is used. The final scenario, CM4 is the one that includes extreme NLOS case. The 
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Table 2.3. IEEE 802.15.3a UWB Channel Model Parameter Definitions 



Parameter 


Definition 


A 

1 L 


A A/PT*?1 Q"P flTTlA/'Pll Tflfp OT pl"N<3TPT*<3 


Ai 


Average arrival rate of multipaths with in clusters (rays) 


F 

1 


JTUVVt;! U-CL-dy IdttUI 1U1 ClLloLclo 


?/; 

Y 


rDW/Pr (iPPflAA TPlpfoT TOT TflA/'C; AA7"lf TllTl f*ln<3TPT*<3 
.L UWC1 U-C-L-Cly LC\i\_,\j\J L 1U1 1 dy o VV 1 ullJ.ll L-lUOLCliJ 


°l 


Standard deviation of the fluctuations of the 
channel coefficient for clusters 


°t 


Standard deviation of the fluctuations of the 
channel coefficients for rays within clusters 


°9 


Standard deviation of the channel amplitude gain 



channel model parameter values for these four channels are given in Table 12.41 

It should be emphasized that since this is a discrete-time impulse response model, 
the time resolution of the system plays a crucial role. Table 12.41 correspond to a time 
resolution of 167 psec which corresponds to a bandwidth of 6 GHz. Therefore, for 
any other time resolution, scaling in time domain is needed for each of the channel 
coefficient realizations. 

Table 2.4. IEEE 802.15.3a UWB Channel Model Parameter Values 



Scenario 


A (1/ns) 


Ai (1/ns) 


T (ns) 


ip (ns) 


a, (dB) 


a c (dB) 


a 9 (dB) 


CM1 


0.0233 


2.5 


7.1 


4.3 


3.3941 


3.3941 


3 


LOS (0-4 m) 
















CM2 


0.4 


0.5 


5.5 


6.7 


3.3941 


3.3941 


3 


NLOS (0-4 m) 
















CM3 


0.0667 


2.1 


14 


7.9 


3.3941 


3.3941 


3 


NLOS (4-10 m) 
















CM4 


0.0667 


2.1 


24 


12 


3.3941 


3.3941 


3 


Extreme NLOS 
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3. FREE-SPACE OPTICAL SYSTEMS 



Free-space optical (FSO) systems are originally created to supplement fiber- 
optical cable systems and their apparatus and techniques originate from fiber-optical 
systems. Digital information in the form of binary signals are sent through building 
roofs or window-mounted infrared laser diode transmitters. Usually, the systems use 
narrow optical pulses to transmit l's and transmits nothing for O's and this type of mod- 
ulation scheme is called On-Off Keying (OOK). Another common modulation technique 
in FSO systems is pulse position modulation (PPM) where there are M pulse slots of 
T p duration. Depending on the position of the pulse, typically equal energy pulses are 
transmitted. In this scheme, synchronization between transmitter and receiver plays 
an important role to determine the starting positions of the pulses. Transmission effi- 
ciency is enhanced by packetizing data which ensures that dividing traffic into packets 
that can be independently sent and received. Also, similar to many fiber-optical com- 
munications, FSO can support wavelength division multiplexing (WDM) which allows 
a single optical path to carry different signal channels ensured that each signal has a 
different wavelength. 

Typically, 850 and 1550 nm laser diodes are the common commercially available 
products in today's technology. In these systems, the information is transmitted by 
modulating the signal intensity denoted by I(t) which is produced in response to an 
input electrical current signal and this process is called intensity modulation. It is 
important to emphasize that the information in these systems are carried in the inten- 
sity of the transmitted signal contrary to the RF systems where either the amplitude, 
the frequency or the phase of the signal or sometimes any combination of these are 
used to transmit information. The generated optical intensity is focused by a lens and 
released out as a beam of light just like flashlights. Although these pulses are focused 
by a lens, the power of the beam still disperses over long distances. At the receiver, 
the transmitted light is focused back onto a photodetector which generates current 
proportional to the intensity incident on its photodetector area and this process is 
called direct-detection. This process inside the photodetector is analogous to square 
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law devices for RF systems where the integral of the amplitude square of the signal is 
taken over a time period. In many cases, the generated current in the photodetector is 
weak and amplification is employed in the receiver end before it is demodulated. By 
this way, optical to electrical conversion of the signals are accomplished. 



Laser 

Diode 

(LD) 



FSO Channel 



Photo- 
detector 
(PD) 



Electrical to Optical Optical to Electrical 
Conversion Conversion 



Figure 3.1. Electrical to optical conversion of signal via laser diode and optical to 
electrical conversion of the signal via photodetector through the FSO link. 

Due to the underlying structure of the channel, types of modulation and detection 
of optical signals are limited to variations on the optical intensity. This imposes the 
non-negativity constraint on the signals to be transmitted and this can be shown as 



I(t) > 0, Vt E 9£ 



(3.1) 



and the above constraint is meaningful in the sense that the transmitted power can 
physically never be negative [31]. Here, it should be emphasized that the received 
signal in the photodetector that is shown by 



y(t) = V I(t)x(t) + n(t) 



(3.2) 



where x(t) is the modulating data and n(t) denotes the noise in the system, can take 
negative values depending on the noise characteristics. 



Another important physical characteristics in FSO systems is beam divergence. 
The transmitted beam diverges by the time it arrives at the receiver. The degree of 
beam spreading depends on the lens diameter and the received energy at the collecting 
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lens decreases with the square of the link distance. Without investigating the atmo- 
spheric conditions, transmit optical power, receiver sensitivity and size of the collecting 
lens, all impose constraints on the range of the communication link for a given data 
rate [32]. 

In order to increase the range of the communication link, the diameter of the 
transmitting lens needs to be increased. By employing larger lens diameters, the power 
incident on the collecting lens is increased resulting in reduced beam spreading. Al- 
though there are physical limitations on how large the lens diameter can be, there is also 
another problem associated with using large lenses. As the diameter of the lens gets 
larger, the beam spreading gets narrower causing the alignment between transmitter 
and receiver to get difficult. Including the building sway effects and thermal expansion 
or contraction of tall buildings, focusing may cause problems due to narrow beams. 
Therefore, many commercial products employ auto-tracking capabilities to both ends 
which brings additional cost and complexity to the system. To deploy tracking capabil- 
ities to FSO systems, movable mechanical platforms or sometimes articulated mirrors 
are used to ensure correct pointing of the transmitted beam. 

Atmospheric weather conditions such as fog, rain and snow are the other factors 
that limit the range of FSO links. Among these weather conditions, FSO systems 
are most susceptible to fog. For example, in moderate dense fog, the optical signal 
loses 90 percent of its strength every 50 meters so for a link distance of 150 meters 
99.9 percent of the transmitted optical energy will be dissipated [52] • Therefore, while 
deploying FSO systems, weather conditions of the deployment cites need to be carefully 
investigated. 

To solve the range versus reliability problem for FSO systems, they are designed 
with limited link lengths to form an optical mesh topology. By this way, when several 
links fail, the spider web like structure may still ensure the signal to be redirected 
through a different path ensuring availability of the service. On the other hand, the 
mesh topology brings complexity to the system. To control the distribution of the 
signals along the paths, network management type systems are needed. 



19 



3.1. Free-Space Optical Channels 

In order to fully design and analyze the properties of a communication channel, 
the propagation characteristics of the path from transmitter to receiver must be taken 
into account. Therefore, electromagnetic propagation of the waves needs to be inves- 
tigated. There are two types of propagation depending on the medium such as guided 
and unguided transmissions. In guided channels, wave guides are used to confine the 
wave propagation from transmitter to receiver. Fiber-optical cables are the primary 
examples for such communication systems. On the other hand, in an unguided channel, 
the transmitter releases the generated field freely into a medium without any attempt 
to control its propagation except to control its antenna gain pattern. This unguided 
channel is often called as space channel and common examples of this channel are free 
space, the atmosphere or the ocean (underwater) [8]. Unfortunately, the characteristics 
of unguided (space) channels primarily depend on the properties of the medium. Free- 
space channels are the simplest type of unguided channels where the medium between 
transmitter and receiver is free-space. 

In FSO links, inhomogeneities induced by the temperature and atmospheric pres- 
sure cause fluctuations in the received light intensity [5] which lead to turbulence in- 
duced fading. Depending on the characteristics transmission path, this type of fading 
may be severe and thus, limit the performance of the communication system by in- 
creasing the link error probability. Also, weather conditions such as fog, rain and snow 
can lead to aerosol scattering and degrade the performance of FSO links and these 
effects are discussed in [33] and [31]. 

In FSO communication systems, there are two parameters that describe turbulence- 
induced fading such as the correlation length of intensity fluctuations, do and the corre- 
lation time of intensity fluctuations denoted by To. These two parameters closely effect 
the link performance by imposing constraints on aperture diameter and bit duration. If 
the diameter of the receiver aperture denoted by Do is larger than the correlation length, 
do then turbulence- induced fading can be reduced significantly by employing multiple 
point receivers in the receiver side and this technique is called aperture averaging [TU] . 
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Since Do > do can not be always satisfied, to mitigate the effects of fading, several 
approaches have been proposed in [5] . These approaches include temporal-domain and 
spatial-domain techniques. To reduce the effects of fading by temporal-domain tech- 
niques is maximum-likelihood (ML) symbol-by-symbol under the assumption that the 
receiver has the knowledge of marginal fading but does not have the temporal fading 
correlation nor the instantaneous fading states. When the receiver knows the joint 
temporal fading distribution but no the instantaneous fading state, the receiver can 
employ MLSD. As spatial-domain techniques, multiple receivers must be employed at 
the receiver to collect the incident intensity. As in radio frequency communication sys- 
tems, receivers need to be separated as far as possible in order to maximize the receive 
diversity gain. Therefore, uncorrelated turbulence induced fading can be achieved. In 
practice, separating receivers may not be made possible due to the large area they 
would occupy and then correlated turbulence fading would occur. Techniques to mit- 
igate the above mentioned turbulence induced fading effects in temporal and spatial 
domain are explained in [5]. 

An important parameter in Kolmogorov theory is refractive index and it is defined 

as [5] 



where n defines the average index and n% denotes the fluctuation component induced 
by spatial variations of both temperature and pressure. The autocorrelation function 
of rii, represents the spatial coherence of the refractive index denoted by r ni (ri, r%). 
its Fourier transform denoted by $> n (k) is called wave number spectrum and it can be 
represented as 



n(r,t)=n + nx{r,t) 



(3.3) 



$ n (k)=F{E[n 1 (n,t 1 )-n 1 (TS,t 1 )]} 



(3.4) 



where F{-} denotes the Fourier transform operation. 



According to Kolmogorov theory, turbulence characteristics can be divided into 
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three regions. 

• Input Range: This region considers the case when the eddy size is greater than 
outer scale of turbulence Lq, ie, (eddy size> L ). In this region, the turbulence is 
anisotropic and the wind shear and temperature gradient introduces energy to the 
turbulence. Since the energy sources may vary, no general formula is applicable 
to this chaotic region [35] . 

• Inertial subrange: This region considers the case when the eddy size is between 
the inner and outer scale of turbulence, ie, (lo < eddy size < L ). In this region, 
the turbulence is fairly isotropic and the kinetic energy of the eddies dominates 
the energy due to the viscosity dissipation. 

• Dissipation range: This region considers the case when the eddy size is smaller 
than the inner scale of turbulence, ie, (lo > eddy size). In this region the energy 
due to the viscosity of eddies dominates the kinetic energy. 

Depending on the ranges, the Kolmogorov wavenumber spectrum can be found dif- 
ferently. Unfortunately, for the input range, the spectrum is still unknown |35j. For 
inertial range, it can be given by 

$ n (K) = 0.033CtftT 11/3 , (3.5) 

where C 2 is the wavenumber spectrum structure parameter and it is an altitude- 
dependent parameter. For dissipation range, the spectrum is considered to be zero, ie, 
<& n (K) = 0. The following formula is often considered to characterize all three regimes 

$ n (K) = 0.033Cl(K 2 + l/L 2 )- 11/6 exp(-K 2 /Kl) (3.6) 

where K m = 5.92// - This equation is referred to as von Karman spectrum. Using 
Hufnagel and Stanley model, the wavenumber spectrum structure parameter can be 
given by [35] 



C 2 n (z) = K z- 1 / 3 exp(-z/z ) 



(3.7) 
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where Kq describes the strength of the turbulence and zq is the effective height of the 
turbulent atmosphere. Several works in literature define the wavenumber spectrum 
structure parameter as 0, EH] 

Cl(z) =0.00594(w/27) 2 (f0" 5 2) 10 exp(2/1000) (3.8) 
+ 2.7 x f0^ 6 exp(-^/1500) + Aexp(-z/1000) 

where z is altitude in meters, v is the root mean square (rms) of the wind speed in 
meters per second and A is the nominal value of C 2 (0) at the ground and its unit 
is m~ 2 / 3 . The wave number spectrum structure parameter, C 2 can vary from 1CT 17 
m - 2 / 3 for weak turbulence regime to 10~ 13 mT 2 / 3 for strong turbulence regime. A 
typical average value is assumed to be 1CT 15 m~ 2 / 3 in [9]. 

This work assumes that the energy of large scale eddies are redistributed without 
loss to eddies of decreasing size until finally dissipated by viscosity [5]. Therefore, the 
size of turbulence eddies may vary from millimeters to meters. Denoting the inner and 
outer scale of the eddies by l and L , respectively, when the link distance L satisfies 
the following condition 

l < \f\L < L (3.9) 

where A is the wavelength of the transmitted signal, 1$ and Lq are defined as above, 
then the correlation length of the fading channel, do can be approximated by [10] 

d w V\L (3.10) 

where a/AL denotes the Fresnel zone of the turbulence. Here, it needs to be mentioned 
that, the above equation does not take into account of aerosol effects such as rain 
and snow. These effects would further degrade the coherence of the optical field and 
lead to decrease in the correlation length of the fading channel. In literature, "frozen 
air" model is used while modelling the behavior of the eddies in the atmosphere. As 
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the name suggests, this model assumes that the eddies present along within the path 
preserve their relative location to each other. 



In FSO channels, depending on the strength of atmospheric turbulence the dis- 
tribution of light intensity is categorized into three regimes such as weak, moderate 
and strong turbulence conditions. All these conditions are represented via distinct 
distributions that accurately model the channel fadings under the associated regime. 
The weak turbulence is represented by log-normal distribution, moderate turbulence is 
described by gamma-gamma distribution and at the extreme case, strong turbulence is 
denoted by negative exponential distribution. These statistical models will be analyzed 
in detail in the following subsections. 

3.2.1. Weak Turbulence Conditions 

The emitted light from the transmitter propagates through large number of el- 
ements of the atmosphere and at times where each element on the propagation path 
causes independent, identically distributed (i.i.d.) scattering and phase delay, the cen- 
tral limit theorem (CLT) can be invoked to represent the marginal distribution of 
log-amplitude fluctuations that will be represented by X and its distribution can be 
shown as 



where E[X] denotes the ensemble average of log-amplitude fluctuations and a\ denotes 
its variance. Commonly, it is assumed that the light intensity represented by the 
random variable, I is related to the log-amplitude X by [5] 



3.2. Statistical Models of Turbulence Regimes 




(3.11) 



I 



7 exp(2X - 2/j, x ). 



(3.12) 
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Then, the marginal distribution of the light intensity can be expressed as 

The variance of log-amplitude fluctuations of plane and spherical waves can be found 
using the following equations [5j |9] 

7/6 r L 

2| 

x I plane 



0.56 (j^J C 2 n {x){L - x) 5/6 dx (3.14) 

spherical = 0.56 (y) ? * £ (%(x) (~) ^ (L - x) 5/6 dx. (3.15) 



where is the wavenumber spectrum structure parameter as described previously, 
and L is the propagation path between the transmitter and receiver. Also, below 
expression is commonly use to relate scintillation index and variance of log-amplitude 
fluctuations such that 

°si = ex P i a x) ~ !■ ( 3 - 16 ) 

Although above expressions can handle varying wavenumber spectrum structure pa- 
rameter, in this work, it is assumed to be constant along the horizontal path. On the 
other hand, in some references such as [H [TTJ, [17] , the log- normal distribution of the 
light intensity has a similar representation which uses the Rytov theory, and it is given 

as 



where aj denotes the aperture-averaged scintillation index. Depending on the wave 
characteristics, the scintillation index for a plane wave and spherical wave can be 
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expressed as 



C" 1 1 plane — © X P 



0.49a? 



® 1 1 spherical exp 



(l + 0.65c/ 2 + 1.11(7 1 12/5 ) 7 6 



0.49a 2 

1 + 0.18rf 2 + 0.56(T 2 12/5 y /6 



+ 



+ 



0.51a 2 (l + 0.69CT} 275 ) 5/6 
1 + 0.9c/ 2 + 0.62c/V 1 12/5 ) 5/6 



0.51a 2 (l + 0.69a 2 12/5 ) 5/6 
1 + 0.9c/ 2 + 0.62c/ 2 cr 2 12/5 ) 5/6 



- 1 



(3.18) 



- 1 



(3.19) 



where a\ = 1.23C 2 fc 7 / 6 L n / 6 and a\ = 0.4a 2 are the Rytov variance for a plane wave 



and spherical wave, respectively, d = ^JkD 2 /(4L), D is aperture diameter, k is the 
wave number given by k = 2ir/X, and the parameters C\ and L are given as above. In 
cases where only a single point detector is used, ie, aperture averaging is not employed, 
d is taken as zero since the propagation path, L, is much larger than the diameter of 
a single detector, d. Finally, although using the relation between light intensity and 
Rytov variance is widely used as in [U [17], it is not considered in this work under 
the weak turbulence conditions but similar expressions invoking Rytov theory under 
moderate turbulence conditions are considered in following section. 



In the following two figures, the close relation between standard deviation of log- 
amplitude fluctuations, o x and path distance, L is shown. Using (I3.14p and A3. 15j) . as 
the link distance varies from hundred meters to a kilometer, the standard deviation 
of log-amplitude fluctuations increase from 10~ 3 to 10 _1 . This general behavior is 
depicted in the figures and the wavelength of the optical transmitted wave is taken as 
1550 nm and C 2 is assumed to be constant throughout the link distance. 
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Figure 3.2. Standard deviation of the log-amplitude fluctuations versus propagation 
distance for a plane wave. The wavelength is taken as A = 1550 nm. 
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Figure 3.3. Standard deviation of the log-amplitude fluctuations versus propagation 
distance for a spherical wave. The wavelength is taken as A = 1550 nm. 
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3.2.2. Moderate Turbulence Conditions 



Under moderate turbulence conditions, both the small and large scale scatter- 
ers are effective and the intensity fluctuations are modelled by the Gamma-Gamma 
distribution described by the following probability density function 



(3.20) 



where T(-) and K n (-) are the Gamma function and the nth-order modified Bessel 
function of the second kind, respectively Small and large scale scattering effects 
that are represented by a and /3, respectively in ( 13 .20 j) are related to the atmospheric 
conditions for a plane wave with aperture-averaged scintillation index by [10] 
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where a\ = 1.23C 2 k 7 / 6 L u ^ is the Rytov variance for a plane wave, 
D is aperture diameter and the parameters C 2 , k and L are given as above. The same 
parameters are related to atmospheric conditions for a spherical wave with aperture 
averaging by 
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(3.22) 



where a\ = 0.4a 2 (ie, a| = 0.5C 2 fc 7//6 L 11 / 6 ) is the Rytov variance for a spherical wave 
and the other parameters, C 2 ,d,k and L are as defined previously. These parameters 
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are related to the scintillation index by ajgj = a 1 + (3 1 + (a/3) 1 and moderate 
turbulence regime is often assumed to be valid for 0.75 < a 2 sl < 1. 

Increase in the link distance L or increase in the wavenumber spectrum structure 
parameter C 2 would result in an increase in the Rytov variance for each type of waves, 
<j\ and a\. Using this fact and the above equations to determine small and large scale 
effects, the following figure is plotted to depict the relation between link distance and 
(a,P,agj) parameter pairs. For small link distances, both small and large scale scat- 
tering effect parameters a and (3 have large values and as the link distance increases 
small scale scattering effects a significantly increases after some point whereas large 
scale scattering effects (3 constantly decrease and approach to unity. For the scintilla- 
tion index, this parameter always increases and reaches a saturation limit as the link 
distance gets larger. This saturation regime is better described by negative exponential 
distribution which is explained in the next subsection. 




2 4 6 8 10 12 14 

Link Distance (km) 



Figure 3.4. Small and large scale scattering effects and scintillation index as a 
function of link distance. The wavelength and wavenumber spectrum structure 
parameter are taken as A = 1550 nm and = 10~ 14 m~ 2 / 3 , respectively. 

3.2.3. Strong Turbulence Conditions 

Under strong turbulence conditions, there are many non-dominating scatterers 
in the medium. Passing from medium to strong turbulence regime corresponds to the 
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case when the small scale scattering effects denoted by a increase but the large scale 
scattering effects approach represented by approach to unity as the link distance L 
increases and under this condition, the light intensity is better modelled as one-sided 
negative exponential distribution which is given by 



cxp 



/ > 



(3.23) 



where I denotes the mean light intensity. This saturation regime corresponds to scintil- 
lation index for the range crf^ > 1. To show that gamma-gamma distribution converges 
to negative exponential distribution in the saturation regime the following figure is plot- 
ted. For the same (a, (3) pairs in the previous subsection, the pdf 's are depicted in Fig. 
13.51 As the link distance increases (3 approaches to unity and a grows unboundedly. 
This behavior is depicted in the figure such that the longest link distance which is 
shown as the red curve has a negative exponential distribution. Aside from this math- 
ematical interpretation, physically in the saturation regime large scale effects are not 
dominant instead there are several non-dominating small scale scattering effects are 
effecting the propagation of light. 



As the link distance increases 




Light Intensity ( I ) 



Figure 3.5. The probability density function of several (a, (3) pairs as a function of 

link distance. 
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4. PROPOSED ULTRA- WIDEBAND SYSTEM AND ITS 

ANALYSIS 



The proposed system model in this work includes two separate sections such as 
free-space optical (FSO) and radio-frequency (RF) system models. Both of the systems 
employ ultra-wideband signalling to achieve very high data rates. The FSO system is 
designed such that it can take binary inputs in both electrical or optical domain. These 
inputs are used to generate optical pulses and transmitted via FSO links over several 
kilometers. The received signals are demodulated in the receiver and they can be used 
right away since they are converted to electrical domain using the photodetectors in 
the receiver. The main advantage of this section of the system is to reduce the time and 
cost to install fiber optical cables for long distances since both the installation time and 
cost to deploy FSO equipment is much shorter and cheaper compared to deploying fiber 
optical equipments. Also, depending on the link design, the received estimates of the 
FSO links can be used to distribute the information inside the buildings. Using simple 
and commercially available UWB transmitter and receiver pairs, the information can 
be delivered to the end-user in the last-step. Due to the severe attenuation of the 
signal in the RF domain, the distance between the transmitter and receiver in the 
last-step is limited to short distance up to 10 meters. For both sections of the system, 
their systems models are presented and their detection error probabilities (DEP) are 
analyzed in the following sections. The FSO system performance is closely related to 
the atmospheric turbulence conditions. Therefore, its performance analysis is divided 
into three parts depending on the strength of the turbulence such as weak, moderate 
and strong conditions. Also, simulations are carried out to verify the derived theoretical 
results. Towards the end of this chapter, channel capacity analysis of the FSO system 
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is presented. 

4.1. System Models 
4.1.1. Free-Space Optical System Model 



The system model shown in Fig. 14.11 is considered and the approaches presented 
in [37] and |38j are adopted. The proposed system employs gain-switched Fabry- 
Perot laser diode (FPLD), a tunable-filter (TF) and an erbium-dope fiber amplifier 
(EDFA) to generate wavelength-tunable optical pulses and is primarily developed for 
UWB-over-fiber links. In the system model, the binary source outputs {uaJ's are time- 
hopping pulse-position modulated (TH-PPM) via Gaussian pulse trains and passed 
through a bias-tee circuit to drive the FPLD into gain-switched operation. The 1550- 
nm FPLD, which operates with a threshold current of 18 mA at 25° C with 0.8 nm 
mode spacing is biased at 16 mA and gain-switched at 4 GHz. The generated optical 
signal is fed to the EDFA which consists of a 980-nm pump laser diode that pumps 
50 mW output power to couple an erbium-doped fiber via 980/1 550-nm wavelength 
division multiplexer (WDM) and an isolator to reduce back reflections and serves both 
as an external-injection source and an amplifier for the FPLD output. The output of 
EDFA is passed through TF which operates in the range from 1527 to 1562 nm. The 
central wavelength of the TF is chosen to be close to that of FPLD output so that the 
system has a single wavelength output before it is sent to FSO channel. The generated 



RF Input Signal 





II,. 


Source 


i\ 

— P(t) J 




/V 

Ui. 





Amp 



DC Bias 



DC , 



FP-LD 



End-User 



RX 



A 
Ul 



UWB 
TX 




Matched Filter 



P «ro\>- P r (T-t) | — CK+K PIN-PD 



Figure 4.1. System model of the proposed communication system over FSO + RF 

links. 
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optical TH-PPM signal can be represented as 

oo N s -1 

x F {t) = ^/E^^ Y p ^ ~~ nTd ~ 3 T f ~ c i Tc ~ dnT v> ( 4 " 1 ) 

n=l j=0 

where E F denotes the pulse energy amplified by EDFA, p(t) denotes the Gaussian 
pulse, {d n } and {cj} are the binary information and pseudo-random code sequences 
for time hopping, respectively. Td, Tf, T c and T p represent the symbol, bit, chip and 
pulse durations, respectively and bits are repeated N s times in a symbol period. 

The FSO channel is described by the unit impulse response h(t) = I(t) + Ib where 
I(t) and h are the instantaneous light intensity and the background radiation whose 
effects are removed at the receiver as in [5], respectively. Also, experiments show that 
the coherence time of the FSO links is large enough to approximate the intensity as a 
constant during the transmission of a frame, i.e., I(t) w /. 

The photo-detector receives the photon flux incident on the detector area and 
produces a current that is proportional to the received photons, providing the optical 
to electrical conversion of the received signal. The conversion coefficient r\ (0 < rj < 1) 
indicates the efficiency of the photodetector. After the removal of the background 
radiation bias rjlfy, the resulting received signal can be represented as 

r F (t)=r)Ix F {t)+n(t) (4.2) 

oo Ns-1 

71=1 j=0 

where n F (t) represents the combined effects of both the thermal noise and the shot 
noise, which can jointly be modelled as an additive white Gaussian noise (AWGN) with 
zero mean and variance No/2, i.e., n F (t) ~ Af(0,N /2). The received signal is passed 
through the matched filter xmf(^) which can be written as 

j oo N 3 -l 

XMFit) = —= Yl p ^ ~ c i T ^ ~ p ^ ~ c i Tc - T p) ( 4 - 3 ) 
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resulting in the sampled output 

f kT d 

y F (kT d ) = / r F (t)x M F(t)dt = ±r ] IyfN s E F + v k . (4.4) 

J(k-l)T d 

Here, the noise term v k = J^^ Td riF(t)xMF(t)dt is also AWGN with v k ~ jV(0, N /2). 
To make decisions on the sampled output values zero-threshold detection is employed 
for the matched filter outputs. The estimates of these binary source outputs, {ut} 's 
are then input to an UWB transmitter to be transmitted to the end-user. 

4.1.2. Radio-Frequency System Model 

In this subsection, radio-frequency section of the system model is considered. It 
should be emphasized that this part of the system constitutes the last-step to the end- 
user. The TH-PPM pulse that is generated by the UWB transmitter can be expressed 
as 

oo N s -1 

x R (t) = ^E~ R Y,J2<l(t-nT d - jT f - c 3 T c - d n T p ) (4.5) 

n=\ j=0 

where E R denotes the energy of the pulse, q(t) denotes the Gaussian monocycle pulse 
train satisfying the UWB regulations. The UWB channel is modelled as 

N K(l) 

h R {t) = Kj2J2 aik6 ( t - Ti - Ti ^ ( 4 - 6 ) 

1=1 k=l 

where k denotes the log-normal fading random variable, iV and K(l) denotes the num- 
ber of observed clusters and multipaths within each cluster, respectively. T; is the delay 
of the I th cluster and aik denotes the channel coefficient of the I th cluster k th multipath 
and it can be expressed as a\ k = pikPik where pik is a Bernoulli random variable taking 
values of ±1 and (3ik is the log- normal distributed channel coefficient. To normalize 
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each channel realization to unity requires that 

N K(l) 

EE#* = 1 - (4-r) 

1=1 k=l 

The amplitude gain k is also assumed to be a log-normal r. v. with the relation 
K = 10 9//20 where g is Gaussian distribution with mean g and variance a 2 g . The mean 
go depends on the average total multipath gain G and expressed as 

lOlnG gjlnlO 
90 ~ InTo 20~ (48) 

where ln(-) is the natural logarithm. G is dependent on average attenuation exponent 
if) by G = Gq/DV, and Go is the reference power gain evaluated at D = 1. Aq = 
10 \og w (ETx / Erx) denoting the path loss at a reference distance D = 1 in dB is 
related to Go by Go = 10 -A °/ 10 . Assuming the transmission of the Gaussian monocycle 
pulse train defined in (14.5p . the received signal can be expressed as 

r R (t) =x R (t)*h R (t)+n(t) (4.9) 

00 N s -1 N K(l) 

= ka/Er E ^ E E P ^ ~ nTd ~ i T f ~ °i Tc ~ dnT P ~ T i~ + n R( t )- 

n=l j=0 1=1 k=l 

Using an UWB receiver, estimates of the signal that passed through FSO and RF links, 
{ufc}'s are obtained. 

4.2. Detection Error Probability Analysis 

In this section, the error probability of the composite system including the in- 
troduced by both FSO and RF links is presented. The total error probability of the 
complete hybrid system is denoted by P e {.lF,lR) an d it can be expressed as 



Pe(lF, 1r) = PeAlF) [1 - PeAlR)} + P^Jr) [1 " ^,f(7f)] (4.10) 
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where P ej i?(7i?), P e .R{lR) an d 7f, 1r are the error probabilities and symbol signal- 
to-noise ratios of the FSO and RF systems, respectively. Since binary TH-PPM is 
employed, the above equation includes errors introduced by both FSO and UWB links 
of the composite system. Although the primary scope of this work is to analyze the FSO 
system performance for different channel regimes but also, the performance analysis of 
the RF channel is briefly investigated and is presented along with the analysis of the 
FSO channel. 

The detection error probability (DEP) for the binary TH-PPM signal model over 
FSO subsystem assumes perfect synchronization between the transmitter and receiver. 
Also, channel state information (CSI) is assumed to be available at the receiver. Thus, 
the detection error probability can be expressed as 



where erfc(-) is the complementary error function and 7^ = Fj[Ep/N ] is the average 
pulse signal-to-noise ratio (SNR) and E[-] denotes the expectation operation. Notice 
that due to the random fluctuations in its amplitude, the light intensity is modelled as a 
random variable whose distribution is dependent on the turbulence region. Expectation 
is taken over the distribution of the light intensity to find the average DEP such that 




(4.11) 




(4.12) 



4.2.1. Error Probability Under Weak Turbulence Conditions 



Under weak turbulence conditions, the average DEP can be computed by taking 
the expectation in (14.121) over the log-normal intensity distribution in (13.131) such as 
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Unfortunately closed form solution of the above equation does not exist, but via Gauss- 
Hermite expansion in [39] which is given by 



/oo n 
g(y)e~ y2 dy « ^2uig(yi) 
-°° i=i 



(4.14) 



where yi and Ui, i = 1, ..,n, are the 2 th root (abscissa) and associated weight, respec- 
tively, (14.13p can be approximated after a change of random variables y = ln(J / To) / \J&cr 
as 
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(4.15) 



As shown in the simulation results, the Gauss-Hermite expansion up to 20 th order 
whose parameters given in [39] provides a close approximation to the actual integral 
in flOD. 



4.2.2. Error Probability Under Moderate Turbulence Conditions 

Under moderate turbulence conditions, the average DEP can be computed by the 
expectation in ( I4.12p over the Gamma-Gamma intensity distribution given in (13.201) 
such as 



^rfc (• /^MM 2 M)^ r^-l 



r»r(/3) 



J^- 1 ^.^ ^-s/aiS/) d/. (4.16) 



Using the formulation in [40] the integrands erfc(-) and K n (-) in (14.161) can be expressed 
in terms of Meijer's G-function and rewritten again as 
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which, using Eq. 07.34.21.0011.01 of [10], results in a closed-form expression such that 
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(4.18) 



The intermediate steps in deriving the above closed-form expression are presented in 
Appendix IB. II 



4.2.3. Error Probability Under Strong Turbulence Conditions 



Using (I4.12p and the distribution given in f)3.23p . the average DEP for strong 
turbulence conditions can be written as 



oo j 

-erfc 



1 — erfc 



(4.19) 



exp 



WVW^fJ V2(r//) 2 iV s7F 
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which again has a closed form expression. 



4.2.4. Error Probability in Coded FSO Links 



The system proposed in Fig. 14. II can also be implemented with a channel encoder- 
decoder pair to protect the information bits transmitted over the FSO channel against 
channel effects. If a convolutional encoder is employed together with the Viterbi de- 
coding on the receiver side, the coded bit error probability p^ ed (jjpj in terms of the 
uncoded error probability is expressed approximately as in pp. 531 of [JT] 

KT M ~ N b [4P e , F ( lF ) (1 - P e , F { lF ))] d ^ 12 (4.20) 

where df ree is the minimum free distance of the convolutional code, and Nb is the 
sum of the Hamming weight of all the input sequences whose associated convolutional 
codeword have a Hamming weight of df ree . 
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4.2.5. Error Probability in RF-UWB Links 



Assuming normalized channel coefficients, i.e., ( 14 .71) holds, the received SNR at 
UWB receiver can be expressed as ^rf = K 2 N s -fR, where 7^ = E[Er/N ] denotes 
the average pulse SNR and results in the following error probability P e ,R{lRi re) = 
|erfc i^sj '7 rk 2 1 '2^j . Averaging out the pdf of P e ,R(lR-> re) over re results in 



P.*(7*) = r-*M\r—) ^exp fj201og lor , o r 



(4.21) 



This can be solved by Gauss-Hermite expansion as in ( 14.141) after a change of variables 
y = (20 log 10 re - <7o)/ ( v 72 ^! ) as 



PeAln) - ^ E ^erfc ( V ™ ) (4.22) 



where yi and Ui are as defined earlier. 



4.3. Simulation Results 



In this section, simulation results for the error performance of the proposed FSO- 
UWB system under weak, moderate and strong turbulence conditions and for the 
hybrid system including the RF section in the last-step are presented. For the simula- 
tions, frame size, repetition rate and the efficiency of the photo-detector are chosen as 
N = 1000 symbols, N s = 2 and 77 = 0.9, respectively. 

For each turbulence condition, it is assumed that C\ is constant through the hor- 
izontal path and for weak turbulence, it is taken as C\ = 5.1 x 10~ 15 . For link distances 
over 2, 2.5, 3 and 4 km, these correspond to (a x , a 2 SI ) = (0.3, 0.095), (0.37, 0.15), (0.44, 0.209) 
and (0.57,0.380), respectively. The analytical approximation derived in (14.151) and the 
Monte-Carlo simulations are shown in the group of first 4 BER curves in Fig. 14.21 
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As seen in the figure, the theoretical derivations, represented with dashed-lines show a 
good correspondence to the simulated performances, depicted in solid lines. 

For moderate turbulence condition, the wavenumber spectrum structure param- 
eter is taken as C\ = 1.76 x 10~ 14 and consider link distances over 2, 2.5 and 3 
km. Using the relation in (13.211) . the (a,/3,cr| 7 ) parameter pairs are (4.16,2.21,0.8), 
(4.00, 1.75,0.96) and (4.05, 1.51, 1.07), respectively. The simulation results in Fig. 03] 
for all three cases are in full accordance with the exact closed form expressions given 
in terms of the Meijer's G-function. 

For strong turbulence, unit mean light intensity is assumed, I — E[I] = 1 and 
a 2 SI = 1. The analytical and the Monte-Carlo simulated results are shown in the figure 
and as expected, the system exhibits the worst performance in this case due to the 
saturation characteristics in this regime. 



Figure 4.2. Comparison of the simulated and theoretical DEP's for weak, moderate 

and strong turbulence fading conditions. 

In Fig. I4.3[ the transmission of the detected optical signals through a LOS RF 
UWB channel in the last-step to the end-user is considered. The parameters for LOS 
channel model, CM1 are taken as a g = 3 dB, A = 47 dB and ip = 1.7 [301 E2] 
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Figure 4.3. FSO + RF UWB System Performance of a link distance of 2 km. For 
weak, moderate and strong FSO turbulence, channel parameters are taken as 
a x = 0.3, (a, /3) = (4.16,2.21) and E[I) = 1, respectively. For RF environment, CM1 

is assumed. 

and the RF link distance is taken as 2 meters. For the FSO link, we assume that 
the signals are transmitted over 2 km at 30 and 45 dB. Under weak and moderate 
conditions, C\ is again chosen as 5.1 x 10~ 15 and 1.76 x 10 -14 which corresponds to 
a x = 0.3 and (a, (3) = (4.16,2.21), respectively. The results are shown in Fig. 14.31 
where the thick dashed curve shows the performance of the RF UWB transmission 
over CM1 channels only. The rest of the curves consider the case where the detected 
bit values are transmitted over the FSO channel. Therefore, these curves represent the 
combined error performance over both the FSO and the RF sections of the transmission 
link. Notice that the hybrid system performance exhibits error floors resulting from the 
errors introduced in FSO links. However, if the SNR over the FSO link is over 30 dB, the 
error floor for weak turbulence is significantly low. This indicates at this or higher SNR, 
under weak turbulence, the hybrid system increases the UWB signalling range to 2 km 
without any sacrifice in the performance. Similar results can be achieved for moderate 
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turbulence conditions, if the FSO channel SNR is 45 dB or more. Unfortunately, for 
the extreme case of strong turbulence, the error floors appear high even at large SNR 
values. 




6 8 10 12 14 16 18 20 22 24 26 

Signal-to-Noise Ratio y R (dB) 

Figure 4.4. (27, 31)g convolutional coded FSO + RF UWB system performance over 

2 km. 

To lower these error floors, we consider the case when a simple convolutional en- 
coder/Viterbi decoder pair within the FSO subsystem is employed. With this addition, 
the errors introduced in the FSO links are reduced to improve the overall system per- 
formance. A rate 1/2 convolutional code with generator (27, 31)§ and Nj, = 2, df ree = 7 
([4~T]. pp. 540) is employed together with a Viterbi decoder. The coded system per- 
formance under each turbulence regime is shown in Fig. 14.41 and only the error floors 
under strong turbulence conditions are shown for comparison purposes with uncoded 
signalling scheme. Notice that the use of even a very simple code reduces the error 
floors significantly. For instance, the error floors for the strong turbulence condition 
at 7^ = 30 dB is 1.5 x 10~ 2 for the uncoded system whereas it is reduced signifi- 
cantly down to 7.4 x 10~ 5 when coding is employed. Table 14.11 shows the error floors 
in both cases for each turbulence condition. Thus, at the cost of additional complexity 
introduced by Viterbi decoder, the error floors of the uncoded system performance is 
significantly lowered. 
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Table 4.1. Error floors of the composite system error performance for the uncoded 



and coded cases. 







Weak 


Moderate 


Strong 


7 = 30 dB 


Uncoded 
Coded 


7.4 x 1(T 7 
9.0 x 10" 20 


1.8 x 10" 3 
6.5 x 10~ 8 


1.3 x 10- 2 

7.4 x 10" 5 


7 = 45 dB 


Uncoded 
Coded 


1.7 x 10~ 13 

4.8 x 10~ 43 


5.0 x 10~ 5 
2.2 x HT 13 


2.5 x 10~ 3 
2.0 x 10~ 7 



4.4. Channel Capacity Analysis 



In his pioneering work in 1948, Claude Shannon introduced the term "channel 
capacity" where the bandwidth and the signal-to-noise ratio together determine the 
quality of a transmission channel [13]. Until that moment, this relation was not seen. 
In his works, the capacity of a noisy channel is defined as the maximum possible 
transmission rate at which reliable communication over the channel is possible [22] . 
For transmission rates that are lower than the channel capacity C, one can always 
design codes that achieve arbitrarily small error probabilities. Hence, for a discrete 
memoryless channel, the average channel capacity can be expressed as 



C = maxI{X;Y) (4.23) 

p(x) 

where I(X; Y) denotes the mutual information between the random variables X and Y 
which are commonly referred as the channel inputs and channel outputs, respectively. 
The probability distribution function denoted by p(x) represents the channel variations. 
The average capacity of a power-limited AWGN channel is given by 

C = ~ log ( 1 + £ ] (4.24) 



2 ° V N 



where P and N denote the power of the signal and noise, respectively and P/N shows 
the signal-to-noise ratio (SNR). In this form, the capacity is in bits/sec/Hz units if the 
logarithm is taken in base 2 and in nats/sec/Hz units if the natural logarithm is used. 
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For a given bandwidth B, the same channel capacity can be represented as 



C = 51og(l + ^) (4.25) 



and in this form, its is in bits/sec units since the bandwidth is used in the equation. For 
different channel variations, the above can be revised to include the effects of channel 
that are represented by h as 

E[C] = J B log (l + h^j p(h)dh (4.26) 

This means that the average capacity of a fading channel is found by taking the ex- 
pectation of the AWGN channel capacity over the fading term. This is reasonable to 
make a comparison between the channel capacity of any fading channel and AWGN 
channel. The former is always greater or equal to the latter capacity and to prove this 
Jensen's inequality can be used such as 

E[C] = j B log [\ + h^j p(h)dh < B log + ^ J hdhj = B log (l + 04.27) 

and using the above equation, it can be concluded that the average capacity of an 
AWGN channel is an upper bound to any fading channel capacity. In the following 
subsections, the average capacity of the proposed free-space optical link is investigated 
under weak and moderate turbulence conditions. 



4.4.1. Channel Capacity Under Weak Turbulence Conditions 



On the propagation path of the transmitted light in FSO channel, there are 
several inhomogeneities in the atmosphere which independently lead to distributed 
phase delay and scattering. Using this and invoking the Central Limit Theorem (CLT) 
then X is Gaussian distributed random variable with mean /i x and variance a\. These 
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inhomogeneities in the atmosphere are related to the received light intensity by 

I = I exp(2X). (4.28) 

Here, it needs to be mentioned that in previous sections the above equation included 
the mean term where in this section it is ignored to include its effects in the closed-form 
expression of the channel capacity. The marginal distribution of the light intensity can 
be represented as 

Taking the expectation over the pdf of received light intensity would result in 

m = r < J.MWzW) dI . (4.30) 
1 ' h 21^1 V \ Sal J ' 

Applying change of variables such that y = In (7/ Iq)—2h x then dJ = Jo exp (y + 2fi x ) dy. 
Hence, 

win r 'o^ ( ,.o. u 



Completing the squares in the exponential term would result in 

/oo 
fv(y)dy 

E[I] = I e 2tM * +2a x. (4.32) 

Consider the case when the fading does not attenuate or amplify the average power, 
then normalizing the fading coefficient such that E[I] = I would require to choose the 
mean of the log-amplitudes to be fix = _ & x [6J. Under this condition, the average 
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capacity of the log-normal channel can be calculated as 

' 1 Jo 21^1 \ J 

The above expresses the capacity in bits/sec units since the logarithm is taken in base 
2 and the bandwidth of the message B is considered in the equation. Applying change 
of variables such that y = (ln(/// ) + 2cr^) / y/Sa^, then (I4.33P becomes 

E [C] = J°° -j= log 2 (l + N slp (ploefV**- 202 *) 2 } exp (-y 2 ) dy. (4.34) 

Integrals of this form can be approximated via Gauss-Hermite expression form as in 
( 14.141) . Then, rewriting (14.341) will result in 

E l C \ = B E 7= lo ^ (l + N s1p (Vo^v^I- 2 ^) j (4.35) 
where yi and Ui are as described in previous section. 

4.4.2. Channel Capacity Under Moderate Turbulence Conditions 

The channel capacity of an FSO link under moderate turbulence conditions is 
calculated by solving 

E[C] = B x r( ^ ( r^ (2) l™ J ^ _1 ln (1 + NsMvI) 2 ) K a ^(2^I)dI. (4.36) 

where the pdf of the light intensity under moderate conditions is used in the integration. 
For the above equation, applying change of random variables such that V = I 2 and 
Meijer's G-function representations of the integrands ln(-) and K n (-), it can be written 
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as 



E[C] = B x 



r(a)r(/5)ln(2) J 



1,1 
1,0 



x ^02 a/3W 



a-P P-a 
2 ' 2 



dV. 



(4.37) 



This form can be further simplified to a closed-from expression such as 



a + P 



E[C] 



B 



a/3 



47ri»rOS)ln(2) \y/rj*N~y^ 



x G 



6.1 
2.6 



(a/3f 



-a-P A- a-P 
4 ' 4 

a-P a-P+2 p-a P~a+2 -a-P -a~P 
4 ' 4 ' 4 ' 4 ' 4 ' 4 



(4.38) 



where the intermediate steps to derive the above closed-form expression are shown in 
Appendix IB. 21 in detail. 



In Figure I4.5[ the average channel capacity of the FSO link is presented. De- 
pending on the turbulence conditions, different fading effects are investigated. Under 
weak and moderate turbulence conditions, the average channel capacity of the FSO 
system are plotted. The wavelength is taken as 1550 nm and the wavenumber spectrum 
structure parameters for weak and moderate turbulence are taken as C\ = 5.1 x 10~ 15 
and C\ = 1.76 x 10~ 14 , respectively. The link distances of 2,2.5 and 3 are considered 
in the figure. As expected, when the system is under weak turbulence condition, the 
average channel capacity is higher compared to the case under moderate conditions. 
This behavior is clearly depicted in the figure. Also, another way to interpret the fig- 
ure is the increase in the link distance results in increase in the scintillation index and 
this decreases the channel capacity. This is due to the fact that scintillation effects 
causing amplitude variations become more dominant to reduce the channel capacity 
as the scintillation index increases. The AWGN channel capacity is also included in 
the figure to show that it is the upper bound and small scintillation index values yield 
close average capacities to this bound. 
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Signal-to-Noise Ratio (dB) 



Figure 



4.5. Average channel capacity under weak and moderate turbulence conditions 

for different link distances. 
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5. CONCLUSIONS 



This thesis outlines the benefits of UWB systems that operate in the unlicensed 
spectrum which can provide very high data speeds and along with the great interest 
from both the research community and the industry, UWB systems clearly have a po- 
tential in various fields of electrical engineering such as WPANs, sensor networks includ- 
ing tracking and ranging applications, vehicular radar systems and imaging systems. 
However, UWB systems that operate in RF channels are limited to short distances. 
To this end, a novel optical UWB system model that uses FSO signals is presented 
and the proposed approach intends to extend the link distance from couple of meters 
to several kilometers. Compared to RF links, FSO systems are not subject to severe 
multipath interference that inhibits the transmission of pulses over long ranges but on 
the other hand, they suffer from atmospheric turbulence fading that severely degrade 
the link performance. 

Along with the necessary channel parameter definitions and their mathematical 
expressions, statistical distributions to model different turbulence conditions in FSO 
channels are also included in this thesis. Using these models, the error probability 
analysis of the proposed FSO system model under weak, moderate and strong turbu- 
lence induced fading conditions are presented. Depending on the system design, an 
extended scenario where the estimates of the FSO link outputs are further distributed 
over RF links to the end-user in the last-step is also considered and its error analysis 
is also included. Both the theoretical and simulation results of two scenarios depict 
that increase in the scintillation index which indicates the strongness of the turbulence 
effects closely influence the system performance. As the scintillation index increases, 
bit error probability of the FSO system significantly increases and saturates for the 
strong turbulence conditions. Also, under weak turbulence regime, the effect of the 
errors introduced in FSO link are almost negligible but as the the turbulence condi- 
tions get worse, high error floors are observed in the system performance. A simple 
convolutional encoder and Viterbi decoder pair in FSO link is proposed as a solution 
to protect the bits and reduce the error floors. When error-control coding is employed, 
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it is observed that the error floors are significantly lowered and the quality of service 
of the system is improved. It should be emphasized that better error correcting codes 
such as LDPC coding can further decrease the error probability at the cost of increased 
complexity. 

Another important system parameter, the average channel capacity of the pro- 
posed FSO system is investigated and its closed form expressions under weak and 
moderate turbulence conditions are also presented in this thesis. Concluding remarks 
of this analysis are that the average channel capacity without any fading is the upper 
bound to the fading channels and increase in the scintillation index gradually decreases 
the FSO channel capacity. Although there are some recent studies on analyzing the 
average channel capacity of FSO links under weak and moderate turbulence condi- 
tions such as in [HUE], the outage channel capacity analysis still remains unanswered. 
Also, the average channel capacity of the composite system consisting of FSO and RF 
links is not investigated yet. Several studies propose hybrid link structure that use 
RF links to back up whenever FSO links are not available. These studies focus on 
narrowband transmission of signals that are not applicable to UWB signalling. The 
link availability drawback of FSO systems arise as a consequence of weather conditions 
on the chosen site. Via employing several FSO transceivers in a mesh topology, link 
availability problem can be reduced. 

The studies on FSO systems will clearly increase the knowledge on atmospheric 
effects and the physics of turbulence regimes, improve better techniques to generate 
optical pulse that advance the technology of FSO devices. It is obvious that FSO system 
performance will continue to improve reaching higher data rates, better accuracy in 
aligning devices and more importantly link availability under all weather conditions. 
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APPENDIX A: Mellin Transform and Applications 



The properties of Mellin transform are very similar to the ones of Laplace and 
Henkel transforms. The Mellin transform of the function f(x) and the inverse Mellin 
transform of F(s) are defined as jH] 

POO 

M{f(x)} = F(s) = / x s ~ 1 f(x)dx (A.la) 
Jo 

i pc+ioo 

M- 1 {F(s)} = f{x)U{x) = ■=—. x- s F(s)ds (A.lb) 

^7T2 J c-ioo 

Some basic properties of Mellin transform are listed in the below table. The following 
Table A.l. Basic Properties of the Mellin Transform 



Linearity 


M{df(x) + C 2 g(x)} = CiF(s) + C 2 G{s) 


Scaling 


M{f(ax)} = (l/a s )F(s), a > 


Translation 


M{x a f{x)} = F{s + a), a > 


Derivative 


M{f^(x)} = (-l)n_g£L F ( S - n), n = 1, 2, 3, . . . 


Convolution 


M{Jo°° f(l)9(u)%} = F(s)G(s) 



table shows some special pairs that are useful while working the Mellin transform. 

Table A. 2. Special Mellin Transform Pairs 



F(s) = M{f(x)} 


f(x) = M-^Fis)} 


it / sin(7rs) 


1/(1 + x) 


T{s)/a s , a > 


g— ax 


r(s) cos(7rs/2)/a s , a > 


cos(ax) 


r(s) sin(7rs/2)/a s , a > 


sin (ax) 


r(«/4) 

8 r(l/2-s/4) 


sin 2 (x 2 ) 
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A.l. Electromagnetic Wave Propagation 

In a random medium, modelling wave propagation leads to difficult integrals that 
cannot be easy solved by classical methods. Even numerical approaches may be hard 
to apply because in most cases, the integrand is defined as the difference between two 
quantities that each may produce a divergent integral or sometimes, the integrand may 
be a product of two functions which one goes to infinity and the other goes to zero. 
For such problems, Mellin transform may provides an easy method. These problems 
commonly arise in wave propagation. The Rytov variance of random fluctuations can 
be found by solving the below integral 



aj = 8n 2 k 2 L I I k$ ti (k) ( 1 - cos ( ) ) d«d£ (A.2) 



i r°° / / t K 2i 

JO V V « 

where k is the optical wave number, L is the propagation path length, £ is the normal- 
ized distance variable and $ n (/t) is the spatial spectrum of refractive- index fluctuations. 
Using Kolmogorov power-law spectrum, $ n (/«) can be expressed as 

$ n ( K ) = 0.033C 2 /sT 11/3 (A.3) 

where C 2 is the refractive-index structure parameter and it is a measure of the strength 
of random fluctuations of the medium. Using the trigonometric identity such that 
2sin 2 (x/2) = 1 — cos (a;) in (1A.2j) and employing (1A.3|) . the below integral represent the 



variance of random fluctuations such as 



aj = lQn 2 k 2 L(0m3)C 2 n J* sin 2 {^f^j ^ (A.4a) 

/ T \ 5/ 6 fl poo 

= 5.211C 2 A; 2 L(^J J £ 5 / 6 d£ J x" 8 / 3 sin 2 (x 2 )dx (A.4b) 



where change of random variables method is applied to the second step. The first 
integral is trivial 

e 5/6 d£ = ^ (A.5) 
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For the second integral, the Mellin transform can be directly applied via Table IA.2I 
such that 



/■oo 

/ x" 8/3 sin 2 (x 2 )dx = M 
Jo 



sin 2 (x 2 );s = -- 



lim 



7t r(s/4) 

I 8~r(l/2- S /4)' 



(A.6a) 
(A.6b) 



Combining above equations, the Rytov variance for plane waves can be represented in 
its final form as 



aj = 1.23C 2 fc 7 / 6 L n / 6 . 



(A.7) 
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APPENDIX B: Meijer's G-Function and Applications 



The Meijer's G-function was defined by the Dutch mathematician Cornelis Simon 
Meijer in 1936 to introduce a general function that includes most of the known special 
functions as particular cases. His definition was able to include many known special 
functions as particular cases. This first definition which was using series expansion was 
made revised by Arthur Erdelyi in 1953 to take integrals in the complex plane. With 
this new current definition, most of the special functions can be expressed in terms of 
the G-function and of the Gamma function. The Meijer's G-function in path integral 
form can be expressed as 



G m ' n z 



ai,-..,a p \ 1 f UT=inbi-3)UU T 0-- a i + 3 ) si mu 

bl ,..., bp a,™ a nUi ra - b 3 + s) nu + i r K- - *) 1 ' } 



where L denotes the integration path, TJ is the product operator for sequences and T(-) 
is the gamma function. The above definition is also referred to as the Mellin transform 
and it holds under the following three assumptions. 

1. m,n,p and q are all integers (ie, m,n,p,q G Z) and < m < q and < n < p, 

2. No pole of any T(bj — s) can coincide with any pole of T(l — — s), which means 
that at — bj j£ 1, 2, 3 for k — 1, 2, . . . , n, 

3. z^O. 

The integral path L has three possible choices are such as 

I. L goes from — ioo to ioo such as all poles of T(bj — s) are on the right of the 
path, while the poles of T(l — at + s) are on the left for j = 1,2, ... , m and 
k — 1, 2, . . . , n. For this path, the integral converges when | argz| < 5n where 



S = m + n- \{p + q) > (B.2) 
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The integral also converges when | argz| = 5n > if 



(q-p) ( o+ g 



) 



> i?e{v} + 1, 



(B.3) 



and a is the real part of integration variable, Re{s} and s approaches to both 
+ioo and — ioo where 



2. L begins and ends at +oo and encircles all poles of T(bj — s) in negative direction 
exactly ones but does not encircle the poles of T(l — a^ + s) for j = 1, 2, . . . , m and 
k = 1, 2, . . . , n. Now, the integral path converges for all values of z, if q > p > 0. 
Another convergent path is when q = p > for \z\ < 1. Also, the path that 
\z\ = 1 converges for Re{v} < — 1 where p = q and v is as defined above in (IB.4[) . 

3. L begins and ends at — oo and encircles all poles of r(l— a^+s) in positive direction 
exactly ones but does not encircle the poles of T(bj — s) for j = 1,2, ... ,m and 
k = 1, 2, . . . , n. Now, the integral path converges for all values of z, if p > q > 0. 
Another convergent path is when q = p > for \z\ > 1. Also, the path that 
\z\ = 1 converges for i?e{v} < — 1 where p = q and v is as defined above in (1B.4j) . 




(B-4) 
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There are several special functions that can be expressed in terms of Meijer's G-function 
such as 



7T 



Z b ~-2 



71 



e z = GJJ I -z 



^cos(2v^) = G 



1,0 

0,2 



sm(2^/z) = G 



1.0 
0.2 



ln(l + z) = G\l I z 



b,b+ h 



Kv{z) — 2^oi2 



1,1 
1,0 







(i 






V —V J 




2 ' 2 / 



(B.5) 
(B.6) 
(B.7) 
(B.8) 
(B.9) 



The following identities are useful to adjust the powers and arguments of Meijer's 
G-function such as 



Z <Jp j( j 



p,q 











(■ 






t) 










(■ 






t 



p,q 



















+ a J 



J £ l+v+(p-q)/2 

(2vr)( fc - 1 ) 5 



G 



km,kn 



kp,kq 1 fck{q-p) 



£1 

k ' 

h. 

k ' ' 



ai + fc — 1 

k ' 



(B.10) 
(B.H) 



a„— k+1 



fcl+fc-1 »£ 



where a in (IB.lOj) is a constant and 5 and t> in (IB.llj) are as defined previously in 
flB.2p and (IB.4p . The following equation is often used while taking the definite integral 
of the product of any two Meijer's G-function. Via convolution theorem, the integral 
of the product of any two Meijer's G-function can be represented in a single Meijer's 
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G-function in its closed form such as 



a p 
t 



1 ^n+fi,m+u / ^ 



dz 



-foi -ft -b 



O4, . . . , Cl n , G? T , Ct n ^i, 



, dp 



-G 



CO 



p+T,q+a 



Olj . . . , CLnj d T) Ct n -|_i, . . . , d 



^1 j ■ ■ • j ^mj , b rn j r \ , . . . , b, 



(B.12) 



Above two equations can be written in a more compact way using ( IB. 101) and ( 1B.121) 

as 



POO 

/ 

^0 











si \ 




Op 1 
















t j 



1 ^ q+a,p+T 



dz 



ci, . . . , c„, 1 - / - h, . . . , 1 - / - b q , c H 



(B.13) 



i, . . . , c CT 



This integral form is quite useful to derive a closed form solution for the error prob- 
ability for the free space optical links under moderate turbulence conditions. In the 
following section, detailed derivation for this condition is presented. 
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B.l. Error Probability Analysis Under Moderate Turbulence 

The integrands erfc(-) and K n (-) are represented in terms of Meijer's G-function, 
the error probability of FSO link under moderate conditions becomes 



PeAlF) - fWW) Jo 2 G °> 2 



a(3I 



a-P 0-a 
2 ' 2 



1 



n 2 $ 

^1,2 



1 

0,| 



dJ 



(B.14) 



where the first Meijer's G-function represents K n (-) and the second one stands for 
erfc(-). In this form, the two terms that are related to intensity, / have different 
powers and it can not be directly used in ( 1B.13j) . Let us represent the term associated 
to erfc(-) similar to the form in (IB. lip as 



G 



2.0 
0.2 



a(3I 



a-P P-a 
2 ' 2 



2tt 



G, 



4,0 
0,1 



2 4 



ol-P a~P+2 P-a P-a+2 
4 ' 4 ' 4 ' 4 



(B.15) 



Using this expression of erfc(-), let us rewrite the integral and also adjust the integration 
variable such that it is same as the one in Meijer term, 



„ , . 2(ap) s ¥ 1 
PeAlF) = -\ L , „ x 



rf[(S+fi+!)_i] ^4,0 

T(a)T(f3) " 8tt 3 / 2 J ' °' 4 



r> 2 fi 

^1,2 



2 1 



a-P a-P+2 p-a P~a+2 



4 ' 4 ' 4 ' 4 



riv s7F (^) 2 


1 " 


2 


0,1 . 



dJ. 



(B.16) 



Let us apply change of variables such that V = I 2 and dl = -^^dV. 



T(a)T((3) 8ttV 2 J 2VV 



n 2 $ 

^1,2 



G, 



1.0 
0,4 



[aP) 2 V 



2 4 



q-fi q-ff+2 g-q ff-q+2 
4 ' 4 ' 4 ' 4 





1 " 


2 


o,| . 



dV. 



(B.17) 
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For the integral in this form, closed form expression in (IB. 131) can be directly applied. 
The power for the square of intensity I 2 is the term represented by / in ( IB. 131) and 
for this integral, it is / = (a + /3)/4, and the roots of the first Meijer's G-function 
integrand are R = {^p, °~f +2 , ^p, " +2 }- The new roots for the result of the above 
integral will be determined by 1 — / — R i: i G {1, 2, 3, 4} and they can be expressed as 



1 _ / _ Rl = i _ 

1 - / - R 2 = 1 - 
1 - / - R 3 = 1 - 
1 - / - R 4 = 1 - 



a + /3 
4 

a + /3 
4 

a + /3 
4 

a + /3 



a — /3\ 2 — a 

4 / ~ 2 
a - /3 + 2\ I- a 



A J 2 
p-a\ 2-/3 



4 J 2 

/3-a + 2\ 1-/3 



(B.18a) 
(B.18b) 
(B.18c) 
(B.18d) 



Also, the arguments of the Meijer's G-function terms in ( IB. 131) are rj* = (a/3) 2 /16 
and = r] 2 N s r YF/2. The argument of the closed form expression will be uj/t]* = 
(8r] 2 Ns'jp) I («/3) 2 • Using the above roots and arguments of Meijer's G-function, (IB. 161) 
can be expressed in its closed form as 



2«+/3-3 



r(a)r(/3)?r 3 /2 



^5,2 



a 2 /3 2 



l-q 2-a 1-/3 2-/3 
2 ' 2 ' 2 ' 2 ' 

0,| 



(B.19) 
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B.2. Channel Capacity Under Moderate Turbulence 



The channel capacity of a FSO link under moderate turbulence conditions is 
determined by solving 



E[C] = Bx 



2(<*0) 



x+0 



I»r(/3) ln(2) J 



J^V- 1 In (1 + N slF ( V I) 2 ) K a _p(2^f3I)dI. (B.20) 



Applying change of random variables such that I 2 = V and dl = ^jdV, then the 
channel capacity can be expressed as 



E[C] = B x 



2(a£) 



A + ^ 



r(a)r(^)ln(2)7o 2W 



T ,£i+£ + I_l 
V 4 +2 1 



In (1 + N.-yrfV) K a _p(2J a{3W)dV. 



(B.21) 



In this form, the Meijer's G-function representation of ln(-) and K n (-) can be introduced 
via flR8l) and ( jR9]) such that 



E[C] = Bx 



r(a)r(/3)ln(2) 



1,1 
1,0 



(B.22) 



x ^02 



q-/3 g-g 
2 ' 2 



dV. 



The arguments of the two Meijer's G-function has different powers of the intensity 
Therefore, let us use ( IB. Ill) to change the argument for K n {-) such that 



a-p P-a 
2 ' 2 



J_ 4,0 / (a/3) 2 F 
27t Cjo ' 4 



16 



Q-/3+2 g-q P-a+2 
2 > 2 ' 2 ' 2 



. (B.23) 
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This form of K n (-) can be used in (IB. 221) to derive a closed form expression such as 



E[C] = Bx 



T(a)T(P) ln(2) 



1,1 
1,0 



X 4tt Go ' 4 



[aftfV 



16 



a-/3 a-P+2 0-a /3-a+2 
2 ' 2 ' 2 ' 2 



(B.24) 



dV. 



For the above integral, previously derived closed form integration formula in (1B.13|) 
can be directly used where its roots are 



(B.25a) 
(B.25b) 
(B.25c) 
(B.25d) 
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-/ 


-w = 
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= 1 


-/ 
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-/ 
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a + ft 
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4 - a 
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-ft 
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—a — 

1 = 

4 


ft 
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a + (3 


i=- a 7 


ft 
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Using these roots, ( IB. 241) can be rewritten in its closed form as 

aft 



E[C\ 



B 



*+/3 
2 



47rl»r(/3) ln(2) \^fN7f F 

(aft) 2 



x G 



6.1 
2.6 



IQNsIfV 2 



t-0 4-a-P 
4 ' 4 



a-P Q-/3+2 @-a f3-a+2 -a-fi -a- 
4 , 4,4, 4,4,4 



(B.26) 
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